Wave turbulence served up on a plate 
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Wave turbulence in a thin elastic plate is experimentally investigated. By using a Fourier trans- 
form profllometry technique, the deformation field of the plate surface is measured simultaneously 
in time and space. This enables us to compute the wavevector-frequency Fourier (k, uS) spectrum of 
the full space-time deformation velocity. In the 3D (k, uj) space, we show that the energy of the mo- 
tion is concentrated on a 2D surface that represents a nonlinear dispersion relation. This nonlinear 
dispersion relation is close to the linear dispersion relation. This validates the usual wavenumber- 
frequency change of variables used in many experimental studies of wave turbulence. The deviation 
from the linear dispersion, which increases with the input power of the forcing, is attributed to weak 
non linear effects. Our technique opens the way for many new extensive quantitative comparisons 
between theory and experiments of wave turbulence. 

PACS numbers: 46.40.-f,62.30.+d,05.45.-a 
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Wave turbulence is a state of waves in non-linear inter- 
action as observed for a large variety of systems includ- 
ing Alfven waves in solar winds [l|, [2[ , ocean waves [3[ , 
non linear optics [4] and superfluids [H[. Similarly to 
the phenomenological theory of hydrodynamic turbu- 
lence, the so-called weak turbulence (WT) theory for 
wave turbulence predicts a Kolmogorov-Zakharov energy 
cascade [6[. This analytical weak turbulence theory as- 
sumes the persistence of the space-time structure of the 
linear waves through the dispersion relation. Very few 
experimental studies have taken place and results show 
only partial agreement with theory 0, H, 0, Fur- 
thermore, almost none of these experiments look beyond 
the analysis of measurement at a single point. Here we 
report the analysis of the turbulence of bending waves 
on a shaken, thin elastic plate, a phenomenon used in 
theatres to simulate the sound of thunder. We are able 
to measure the fully resolved space-time dynamics of the 
deformation of the plate and we show that the energy is 
localized on a line in the wavenumber-frequency plane of 
the Fourier spectrum. This confirms the persistence of 
the space-time structure of waves which is the premise 
of weak turbulence theory. In addition, our system dis- 
plays the phenomenology described by the theory and yet 
some of its predictions are not quantitatively fulfilled: the 
non-linear shift to the dispersion relation and the power 
spectrum do not obey the predicted scaling laws. 



ration in the time evolution of the wave amplitude com- 
pared to the wave period and it provides a natural clo- 
sure of the hierarchy of cumulants derived from the wave 
equation @, [ll|. In contrast, no such closure can be 
exhibited for hydrodynamic turbulence. In particular, 
the WT theory of wave turbulence leads to a kinetic 
equation for the evolution of the energy spectrum of the 
waves. Stationary solutions are exhibited which corre- 
sponds to the Rayleigh- Jeans spectrum for systems in 
equilibrium and the Kolmogorov-Zakharov energy cas- 
cade for non-equilibrium systems. This prediction of the 
power spectrum density of the wave amplitude has been 
derived in many cases such as non linear optics, super- 
fluids, gravity-capillary water waves, sound waves, Alfven 
waves, plasmas, oceanography, semiconductor lasers and 
bending waves in elastic plates @, [ill, EH • There are a 
paucity of experiments specifically designed for wave tur- 
bulence and of those, most concern surface waves on liq- 
uids [7|, [8[ • Our wave system consists of a thin steel plate 
on which elastic bending waves are excited by an electro- 
magnetic vibrator. The dynamics of the plate follow the 
Foppl-Von Karman equations for the deformation: 
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The theory of WT relies on the assumption of weak 
non linearity of waves. The latter induces a scale sepa- 



where p is the density, ( the plate deformation, h the 
plate thickness, E the Young's modulus, a the Poisson 
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ratio, A the Laplacian operator, x the stress function and 
{.,.} is a bilinear differential operator [12]. Linearizing 
the first equation in (pQ) provides the linear dispersion 
relation 
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The non linear term in ([T]) is due to the stretching of 
the plate and it is cubic in the wave amplitude. The 
WT theory has recently been applied to this case [T^ 
and predicts a space Fourier spectrum of the amplitude 
of the waves 
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where P is the average power input in the system from 
the applied forcing, C is a number and fc* is a cut-off 
frequency. The one-point spectrum of the waves has been 
investigated experimentally 0, [Io[ and has been shown 
not to obey the WT prediction in particular in its scaling 
in P : 
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Nevertheless, it displays a turbulent-like behaviour , 
i.e. a broadband spectrum, and the question is raised 
whether or not the disagreement with the theory is due 
to strongly non linear structures, to boundary condition 
effects or to some dissipative mechanism. 

A sketch of the experimental setup is shown in fig.[]Ja). 
The plate is made of stainless steel and its size is 2 m by 
1 m and 0.4 mm thick. Its is bolted on a I beam by 
one short end and is hanging under its own weight. An 
electromagnetic vibrator is anchored 40 cm from the bot- 
tom of the plate and excites the waves at 30 Hz with a 
varying amplitude. The vibrator is fitted with a FGP 
sensors force probe and Briiel & Kjaer accelerometer to 
measure the input power P. A Fourier transform pro- 
filometry technique TH, 3] gives access to the temporal 
evolution of the deformation of the plate measured over 
a significant portion of its area. The principle is the fol- 
lowing: a sine intensity pattern I(x,y) oc sin(27n//p) is 
projected on the surface of the plate by a videoprojec- 
tor. The pattern is then recorded by a Phantom v9 high 
speed camera. The deformation of the plate induces a 
phase shift of the pattern recorded by the camera. The 
deformation of the plate is recovered by a 2D phase de- 
modulation of each image in the movie [TH, [HI- Movies 
are recorded either with 1000 2 (resp. 800 2 ) pixels at 1300 
(resp. 2600) frames per seconds (fps). The configuration 
and the processing is similar to that of Cobelli et al. 
with a distance of L — 193 cm from the projector to the 
plate and a distance of D = 35 cm between the opti- 
cal axes. The normal velocity of the plate is obtained 
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FIG. 1: (a) Sketch of the wave turbulence experiment. The 
specifically designed 2D mechanical system is made of a 
2 x 1 m 2 , 0.4 mm thick stainless steel plate held vertically 
and set in motion by an electromagnetic vibrator at 30 Hz. 
The Fourier Transform profllometry is based on the projection 
of a sine intensity pattern by a high definition video projec- 
tor. The deformed image is then recorded by a high speed 
camera (1300 or 2600 fps). (b) Example of measurement of 
the deformation velocity on a 63 cm by 62 cm area. 



by differentiating the deformation movie in time. The 
field of view is about 71 2 cm 2 at 1300 fps and 62 2 cm 2 
at 2600 fps. The spectra are calculated by performing 
a multidimensional Fourier transform without applying 
any windowing to preserve the localization of the energy 
in the Fourier space. An example of the normal velocity 
of the plate is displayed in fig. [TJ 

The full space-time Fourier spectrum (shown in fig. 
|2[a)) of the deformation i£(k, uS) (a function of both the 
wave vector k and the frequency uo) is constructed from 
the movie of the deformation velocity. The striking fea- 
ture is the localization of the energy in the vicinity of a 
surface showing that the motion is a non linear superpo- 
sition of waves following a dispersion relation uo = /(k), 
close to the linear dispersion relation. This is the first 
experimental observation of such a space-time spectrum 
in wave turbulence. In addition to the full space-time 
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FIG. 2: (a) Space-time spectrum E(k, ui) of the deformation 
velocity (colors are log scaled). The cuts are located at k x = 0, 
k y = and co = 5000 rad/s. The energy is localized on a 
surface in the (k, ui) space which confirms that the turbulent 
motion is due to a non linear superposition of waves following 
a dispersion relation, (b) Spectrum E(k, ui) obtained from 
E(k,cu) by integrating over the direction of k. Insert: space 
spectrum E(k) computed from E(\t,u) by summing over the 
frequencies. Contours are log scaled in both plots. 



spectrum, we can analyse the space spectrum E(k), as 
displayed in fig. [2][b). The isocontours for large wave 
numbers are circles, revealing the isotropy of the spec- 
trum in this regime. The anisotropic response to the 
forcing is visible at low wave numbers. This behaviour is 
expected in the phenomenology of the Kolmogorov cas- 
cade of energy and is evidenced here: After a few steps 
in the cascade, the anisotropy of the forcing is forgotten 
down to the small scales at which the dissipation is dom- 
inant. Owing to the isotropy of the spectrum, in fig.|2][b) 
we show the spectrum E(k = ||k||,o;) obtained by inte- 
grating E(\a, ui) over all the directions of the wave vector. 
The localisation of the energy appears as a line in the 
(fc, ui) plane. The width of this line is close to the inverse 
of the image size, which indicates that the localization of 



the energy in our measurement is actually limited by the 
resolution of the Fourier transform due to the finite size 
of the plate. At low frequency, the injection of energy 
corresponds to a peak on the energy line: even though 
the forcing is localized in space, its monochromatic na- 
ture (at 188 rad/s) makes it local in the (fc, ui) plane. The 
forcing operates effectively at low frequency and at low 
wavenumber as is expected in the phenomenology of the 
Kolmogorov- Zakharov cascade. 
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FIG. 3: Non linear dispersion relation u){k) computed from 
the line of maximum energy in the space-time spectrum of 
the plate deformation velocity for various input power P of 
the forcing (from bottom to top P 1 / 2 = 1, 2, 3, 4 in arbitrary 
units) . The dashed line is the linear dispersion relation J3]) . A 
systematic shift is observed which increases with the forcing 
power. The top insert shows the deviation from the linear 
dispersion. At high wave numbers, the shift is seen to be in- 
dependent of k. The bottom insert shows the shift normalized 
by P 1/2 . 



The concentrated line of energy in the spectrum 
E(k,ui) allows the dispersion relation to be computed; 
it is extracted by computing the position of the crest of 
the energy line at each frequency and is displayed in fig. [3] 
for various values of the forcing P. The dispersion rela- 
tion remains close to the linear dispersion relation with a 
small but systematic shift. This provides strong evidence 
that our system is indeed weakly non linear. Thus, the 
quantitative disagreement between the one point spec- 
trum and the WT theory prediction 0, cannot be at- 
tributed to the existence of strongly non linear structures. 
Instead, the disagreement is proposed to be attributed 
to the dissipative mechanisms which are believed to ex- 
ist at all scales rather than being present at only small 
scales; hence the Kolmogorov- Zakharov cascade is "leak- 
ing" £3. 

Figure [3] allows us to quantify the departure of the ob- 
served dispersion relation from the linear one. Notably, it 



4 



is shown that the correction increases with the power in- 
put P with a behaviour close to P 1 / 2 behaviour. It is also 
observed to be constant at high wavenumbers so that the 
various dispersion relations are parallel to one another. 
This behaviour is of particular interest when compared 
with the WT prediction. Indeed, it is expected that the 
departure of the dispersion relation from the linear one 
has a power law scaling in P which is identical to the 
scaling for the energy spectrum [16]. Our experimen- 
tal measurements confirm this prediction: the exponent 
close to 1/2 in P is common to both departure from lin- 
earity in the dispersion relation and also to the energy 
spectrum. 
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FIG. 4: Test of the validity of the change of variable k <-> uo 
using the dispersion relation. Solid blue line: direct mea- 
surement of the time spectrum E(uj) - solid red (resp black 
dashed) line: E(co) computed from the space spectrum E(k) 
using the nonlinear (resp. linear) dispersion relation to change 
variable. Insert: same for the space spectrum E(k). Solid blue 
line: direct estimation of E(k) - red (resp black dashed) line: 
E{k) computed from E(uj) by the change of variable using the 
nonlinear (resp. linear) dispersion relation. 



The application of the WT turbulence theory is often 
restricted to the prediction of the Kolmogorov-Zakharov 
space spectrum E(k). Although some experiments di- 
rectly measure the space spectrum [l?], [H|, ^ * s often 
easier to measure the motion at one given point as a 
function of time. In that case, only the time spectrum 
E(uo) can be estimated. To compare with the theory, the 
space spectrum is determined by using the dispersion re- 
lation to obtain E(k) from E(uS) (JfiEl- This approach 
was used to deduce the law in Eq. ([5]) [1, [To[ . We can in- 
dependently estimate both E(ou) and E(k) directly. We 
can assess the validity of the change of variables tech- 
nique above (via the linear or non linear dispersion rela- 
tion). The comparison of the various cases is shown in 
fig. [U Both dispersion relations allow us to reproduce 
the inert ial range, with a better agreement shown when 



using the non linear relation. The large time or length 
scales are well reproduced only when using the non lin- 
ear dispersion relation. This validates the usual change 
of variables when the non linearity is weak. 

Our experimental approach of wave turbulence reveals 
the main features of the weakly coupled waves that can 
be usefully compared with weak turbulence theory. Un- 
til now, only the spectrum E(k) has been compared with 
theoretical prediction. However, WT theory can go far 
beyond spectra predictions: it gives quantitative predic- 
tions for multipoint statistics. Our present study con- 
firms and quantifies the weakly nonlinear behaviour of 
the waves comprising the turbulent cascade. It confirms 
that the scaling law in the supplied power P is the same 
for the departure from the linear dispersion relation and 
also for the energy spectrum. Overcoming the discrep- 
ancy between experiments and theory claimed previously, 
we have shown some agreement between experimental re- 
sults and WT theory. We anticipate that this experiment 
will allow precise and quantitative comparisons with the- 
oretical investigations of wave turbulence @, 0, Hi] of 



prime importance for the large number of turbulent sys- 
tems in which extensive measurements are out of reach. 
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